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We study the scattering of solitons in the nonlinear Schrodinger equation on local inhomogeneities 
which may give rise to resonant transmission and reflection. In both cases, we derive resonance 
conditions for the soliton's velocity. The analytical predictions are tested numerically in regimes 
characterized by various time scales. Special attention is paid to intermode interactions and their 
effect on coherence, decoherence and dephasing of plane-wave modes which build up the soliton. 
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The scattering of solitons by various scatter- 
ing centers in the nonlinear Schrodinger equa- 
tion leads to resonant transmission and reflec- 
tion if the soliton velocity matches certain res- 
onance conditions. By assuming that the soliton 
is composed of a weighted superposition of modes 
(i.e. a wave packet) different scattering regimes 
are observed depending on the ratio of the dura- 
tion of the scattering event and the characteris- 
tic mode-mode interaction time due to nonlinear- 
ity. Resonant transmission does not suffer from 
mode dephasing, while resonant reflection (Fano 
resonance) does. Consequently transmission res- 
onances are observed independent of the scatter- 
ing regime, while the opposite holds for Fano re- 
flection resonances. 



INTRODUCTION 

Transmission of linear waves through inhomogeneous 
media is a topic of wide interest 0> Ei 01 & @ ■ m the 
one-dimensional case, two basic resonant features of wave 
scattering are known. One of them corresponds to the 
wave travelling above a barrier or a potential well. Prox- 
imity of the wave's frequency to that of a standing wave 
captured inside the scattering potential leads to a reso- 
nance, i.e., the transmission may be strongly enhanced in 
a vicinity of the resonance 6]. In fact, the transmission at 
resonance isperfect, which is used in Fabry-Perrot inter- 
ferometers |j| . The connection of the resonant scattering 
to the Levinson's theorem and the existence of bound 
states in potential wells is discussed, e.g., in Ref. @- 

Resonant reflection is possible too, as a consequence of 
the Fano resonance and it occurs whenever a wave 
may choose between two scattering paths, which finally 
merge into one exit. Destructive interference may then 
lead to total reflection. An adequate explanation goes 



through the coupling of a local Fano state to a contin- 
uum 0, 0, El • The wave may then again either pass 
directly through the continuum or by visiting the Fano 
state. The resonance condition is provided directly by the 
matching of the wave's frequency to the eigenfrequency of 
the Fano state. More than one wave paths are typically 
also generated by time-dependent scattering potentials. 

Fano resonances imply the coherence of the wave 
phases throughout the scattering process, hence any sig- 
nificant dephasing effects will suppress the resonance. 
At the same time, the resonant transmission mechanism 
does not rely on phase coherence, therefore dephasing 
does not destroy it. 

Similar situations can also be observed in nonlinear 
systems which possess continuous or discrete transla- 
tional invariance. For instance, nonlinear lattice mod- 
els support time-periodic spatially localized states in the 
form of discrete breathers [EJ. If small-amplitude plane 
waves are sent towards such a breather, it acts as a time- 
periodic scattering potential. The temporal periodicity 
of the potential leads to excitation of many new scat- 
tering channels, shifted by multiples of the breather's 
fundamental frequency relative to the frequency of the 
incoming wave 13]. Even if all the additional channels 
are closed ones, i.e., they do not match the plane- wave's 
frequency, they may generate localized Fano states which 
resonate with the open channel and thus lead to Fano- 
resonant backscattering [l^j . In addition, resonant trans- 
mission can also be observed in this case [l5|. 

In this paper we study the transmission properties of 
small-amplitude solitons (rather than plane linear waves) 
in the discrete nonlinear Schrodinger (DNLS) equation 
with two types of scattering centers. The first one is 
a Fano-defect center, which is an extra level coupled to 
the DNLS equation at one site of the underlying lattice. 
This Fano center may actually be a strongly localized 
discrete breather of the DNLS model ^ij. The second 
type is a two-site impurity which gives rise to two bound 
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states. Our goal is to predict resonant backscattering 
and resonant transmission of the small-amplitude soliton 
impinging onto these defects. To this end, we consider 
the soliton as a superposition of plane waves or modes, 
while the nonlinearity leads to mode-mode interactions, 
which may or may not cause dephasing effects in the 
course of the scattering process. We will find values of 
the velocity of the incident soliton at which the resonant 
transmission and reflection are possible. 

The predicted effects can be observed in any physi- 
cal system which is modelled by the DNLS equation, the 
most straightforward ones being arrays of nonlinear op- 
tical waveguides. These may be realized as set of parallel 
cores fabricated on a common substrate , or a virtual 
array induced in a photorefractive material illuminated 
by a set of parallel laser beams |l7j . 



RESONANCES IN THE SCATTERING OF 
LINEAR WAVES 

Before considering the scattering of solitons, we will 
set the stage, using two simple models which make it 
possible to observe resonant transmission and reflection 
in linear wave scattering. In both cases, we will use the 
linear discrete Schrodinger equation as the unperturbed 
system. Later on, we will add nonlinearities to make 
propagation of solitons possible. 

Resonant transmission 

Resonant transmission due to bound states is possible 
if the scattering potential supports at least one bound 
state. To be flexible, we take a lattice model which allows 
for two or one bound states, depending on its parameters: 

l(j) n = C((f> n -x + <pn+l) + Sn.O^n + 8 nt l€i<f) n . (1) 

Here 4> n is the complex scalar dynamical variable at the 
n-th site, a real constant C controls the strength of the 
inter-site coupling, and two diagonal defects are set at 
sites n = and n — 1 with strengths eo an d ei, respec- 
tively. In the absence of the defects, Eq. Q has exact 
plane-wave solutions, <f> n = exp [i(ui q t — qn)], whose fre- 
quency satisfies the dispersion relation 

uj q = — 2C cosq . (2) 

To find the transmission coefficient T in the presence 
of the two diagonal defects, we impose proper boundary 
conditions and obtain, after some algebra (see, e.g., Ref. 
flif'). for the case of e = t\ = e, 



e 2 (e - 2C cos q) 2 + 4C 4 sin 2 q ' 



Perfect transmission (T — 1) is possible in the case e < 
2C under the resonance condition (at the special value 
of q), 

q' = cos- 1 (e/2C*) . (4) 

For small values of the impurity strength, e <C 2C, the 
perfect transmission takes place around q — it/2. With 
e approaching 2C, the value of q' moves towards zero 
(see FigEJ. At e = 2C, the perfect transmission occurs 
exactly at the edge point of the spectrum, q' = 0, and 
disappears for larger values of e. Therefore, it is possible 
to control the location of the resonance by varying the 
impurity strength e. 




FIG. 1: The transmission coefficient T versus q for the case 
of the resonant transmission, as per Eq. with C = 10 and 
e = 19. 



Resonant reflection 

In order to produce a Fano resonance, we consider the 
Fano- Anderson model, which amounts to adding an ad- 
ditional local Fano degree of freedom <p, with the eigen- 
frequency (energy) E, to the linear discrete Schrodinger 
equation. The dynamical variable ip is coupled to the 
lattice field <f> n at the site n = 0: 

i(j) n = C(0„_i + 4>n+l) + €ip8 n a , 

i<j> = —Etp + £0o ■ (5) 

When e = 0, the system decouples into the free wave with 
the spectrum uj q = —2C cos q and an additional localized 
level with the energy E. The value of E is chosen so that 
it lies inside the continuous spectrum, i.e., \E\ < 2C. 

If e ^ 0, the Fano defect interacts with the continu- 
ous spectrum locally. To solve the linear system iJJj}, we 
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substitute 

4> n = A n exp (iujqt) , <p = B exp (iu> q t) 
and obtain 



OJqA n 
-LU q B 



C{A n -x+A n+x ) 
-EB + eA . 



eBS n o , 



(6) 



We use the second equation in (JSJ to eliminate B in favor 
of Aq, then the first equation yields 



- UJ q A r , 



C{A n ^+A n+1 )-e 2 (uj q -E) ^o^o- (7) 



Equation (JJJ amounts to the presence of a resonant scat- 
tering potential, which depends on the frequency of the 
incident wave. Moreover, if the energy of the addi- 
tional level is located inside the continuous spectrum, 
| i^l < 2C, as it was assumed above, there is a value 
qp = cos -1 (-E/2C) at which the denominator of the 
last term in Eq. (0 vanishes, according to Eq. (J2J). At 
this value of q, complete reflection, T = 0, takes place, 
as the scattering potential becomes infinitely large. 

After some algebra (see, e.g., Ref. |l8|), the transmis- 
sion coefficient T for the linear model based on Eqs. (JSJ 
can be written as 



T = 



1 



-, -l 



4 sin 2 qC 2 (E~ujq) 2 



(8) 



A typical dependence of T versus the incident wavenum- 
bcr q is shown in Fig. EJfor E = 0. In the case of E = 0, 
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FIG. 2: The transmission coefficient as given by Eq. ©, 
versus the wavenumber q of the inident wave. The parameters 
are E = 0, C = 10 (i.e., qF = tt/2), and e = 4. 

the transmission coefficient T{q) vanishes at qF = tt/2. 
The additional vanishing of T at the band edges q = 
and q = tt is related to the fact that the group velocity 
of plane waves is equal to zero at these points, and will 
not be discussed below. 

If e 2 / (2C 2 ) < sinq, then T w 1 at all values of q 
except for close to q — qF, as T(qp) = 0. The width of 



the resonance in the q-space is defined by the distance 
between points at which T — 1/2, or, as it follows from 
Eq. ©, 



4C 2 (E -toq) 2 sin 2 q 
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By expanding uj q around qF, 
substituting this into Eq. ©, 

5q- 



u q ~ LUq F 

we find 



(9) 

u>'q(qF)Sq, and 



(10) 



4C 2 sin 2 g F 

The width of the resonance is twice the expression JTUJ, 



Aq 



2C* 2 sin" q F 



(11) 



The above considerations are valid also for the model 
based on the continuous linear Schrodinger equation 
augmented with a corresponding 5-like scattering term. 
Quite a different situation takes place if the local scat- 
tering term, added to the linear Schrodinger equation, 
is nonlinear - in that case, the solution to the scatter- 
ing problem is not unique, and may be subject to local 
modulational instability |l9j. 

SOLITON SCATTERING 

In order to consider a possibility of resonant trans- 
mission or reflection in the scattering of a soliton on a 
local defect, we first add nonlinear terms to the linear 
Schrodinger equation [see the first equation in the sys- 
tem JSJ)], thus arriving at the DNLS equation, 



i(f) n = C((f> n -1 + (f>n+l) + \\4>n\ 2 <j>n 



(12) 



Small-amplitude (hence, broad) moving solitons in Eq. 
()12J) are approximated well by the corresponding solution 
to the continuous NLS equation, 



sech 



2C* 
77A 
AC 



exp 



./ V 
l{ 2C n 



2)Ct) 



(n - Vt) 



(13) 



where V is the velocity of the soliton, lo s = 
(4C 2 ) (V 2 — r) 2 A 2 /A) is its intrinsic frequency, and 77 
is the amplitude [20| . 

The velocity V determines a central wavenumber of 
the soliton's Fourier transform, 



q c = VI (2C) 



(14) 



The soliton may be considered as a superposition of lin- 
ear plane waves with wavenumbers taking values in an 
interval around q c , the width of this interval being 



Aq c = In 



(2 + V5) 



(15) 
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This spectral width is to be compared with the width of 
the corresponding resonances [e.g., with the expression 

CD]- 

There are two characteristic time scales, which describe 
the scattering of the soliton by the resonant defect. 
One of them is the time of the soliton-defect interaction, 
which is estimated as the soliton's width 1/Aq c in the 
coordinate space, divided by its velocity: 



r int = (Aq c V)~ 



(16) 



The interaction between waves composing the free soliton 
defines the second time scale. It can be estimated as 
the time of dispersion of the wave packet 1|13[) in the 
linearized equation (|12|l . with A = 0, which yields, similar 
to Eq. C2J> a result 



(Aq c Av g 



(17) 



where V g (q) 
and 



2C sin q is a group velocity of the waves, 



Av g = Vg(q c +Aq c /2)-v g (q s -Aq c /2) = 4 sin (Aq c /2) cos q c 

(18) 

is the relative velocity between faster and slower waves 
composing the soliton. 

The interaction between the plane-waves constituents 
does not play a significant role during the scattering of 
the soliton if 



Tint Tdisp : 



(19) 



hence, under this condition, the soliton may be consid- 
ered as a set of noninteracting plane waves while it suffers 
scattering on the defect. The transmission of each wave 
component is then determined by the corresponding coef- 
ficient for the linear model, see the previous section. If, in 
this regime, Aq c is sufficiently small in comparison with 
the width of the transmission or reflection resonance for 
the plane waves, then all the waves composing the soliton 
will be resonantly transmitted or reflected, provided that 
velocity V matches the resonance condition. 

When Ti n t S> Tdisp, the wave- wave interactions become 
important during the scattering process. Since these in- 
teractions may lead to dephasing of individual waves, 
we expect that the resonant transmission will not be af- 
fected, while the resonant Fano reflection will be sup- 
pressed, as it relies on keeping the wave phase coherence 
in the course of the scattering process. 



We performed numerical simulations of Eqs. (|20|l using 
the fourth-order Runge-Kutta method. The total num- 
ber of sites was large, N = 2000. Time was measured in 
dimensionless units. 

The soliton is launched from the left with a positive ve- 
locity V. After the interaction with the Fano defect, the 
soliton is typically found to be split into two soliton-likc 
fragments which move in opposite directions, see Fig|3] 
The transmission coefficient can be found from numeri- 
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FIG. 3: Profiles of the wave field in the initial state (a) and 
after the interaction (b). The Fano defect (marked by the 
shaded bar) is placed at n = 0, and the direction of motion of 
the pulses is indicated by arrows. The parameters are: E = 
(i.e., q F = f ), C = 10.0, and e = 4.0, 77 = A = 1.0. After the 
interaction, two soliton-like pulses are observed, which move 
in the opposite directions with equal absolute values of their 
velocities. 



cal data, using the conservation of the norm E |</> n | 2 in 

n 

the DNLS equation (another conserved quantity is the 
Hamiltonian) . The transmission coefficient is then de- 
fined as the ratio of the norms of the transmitted wave 
packet and initial soliton: 



*\|2 



T = 



E \Mt 

n>0 

EI^(o)l 2 



(21) 



Resonant reflection of solitons 

A nonlinear model which may give rise to the resonant 
reflection of the soliton is based on a straightforward gen- 
eralization of Eqs. 

i(j) n = C((j) n -i + </> n +l) + \\4>n\ 2 (t)n + t^SnO , 

iip = —E(p + €(j) . (20) 



where we choose t* » T in t- 

In the case shown in Fig. |3 we chose qp = tt/2. In this 
case, the dispersion time l(T7)l diverges, as it follows from 
Eq. (|18[) . and in the first approximation the soliton does 
not disperse at all, even in the linear system. Further, 
the condition Aq c < Aq, which implies that the soliton's 
spectral size is smaller than the width of the resonance 
[see Eq. CD ], hence all the spectral components are 
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expected to be resonantly reflected, is 



To satisfy this condition, other parameters were taken as 
7] = X = 1, C = 10, and e = 4. 

The numerically computed transmission coefficient is 
shown, as a function of the soliton's velocity, in FigQ] 
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FIG. 4: The transmission coefficient for the scattering of the 
soliton on the Fano defect, defined as in Eq. 1211 and found 
from direct simulations, versus the soliton's velocity V. Pa- 
rameters are the same as in Fig. |H] The minimum is located 
at V F = 2Cq F - 

There is a minimum in the dependence of the trans- 
mission on soliton velocity exactly at Vf = 2Cqf, as 
it is predicted by Eq. ltT4*|) and the concept proposed 
above, according to which the soliton may be treated as 
a spectrally narrow packet of quasi-linear waves, hence 
the resonant backscattering should take place when the 
central wavenumber q c coincides with the linearly pre- 
dicted Fano-resonance wavenumber, qp . The minimum 
does not reach zero, which is explained by the finite spec- 
tral width Aq c of the soliton, i.e., wave components which 
build up the soliton do not fully backscatter. It should 
be stressed that the close proximity of the numerically 
found minimum to the anticipated point predicates upon 
the conditions l|22[) and l|19f) . but it is not related to the 
fact that the particular example displayed in Fig. 0]has 
qF = 7r/2. In fact, numerical results are very similar at 
other values of qF- 

To study the case when the wave- wave interactions are 
important during the time of the interaction between 
the soliton and the defect, i.e., Ti nt >Tdi S p, the disper- 
sion time was changed by shifting the soliton's central 
wavenumber q c to the spectrum's edge, keeping all other 
parameters fixed. The Fano resonance is observed in the 
simulations until Tint ~ Tdi sp - With the further decrease 



of ri n t, the Fano minimum in the T(V) curve becomes 
less pronounced, and disappears at some critical value 
of q c . For instance, we have found numerically that for 
C = rj = X = e = l, when the condition of the small- 
ness of rdi S p is fulfilled, the Fano resonance disappears at 
g c »1.0. 



l 



f— 




FIG. 5: The dependence of the soliton's transmission coeffi- 
cient in the case of the scattering on the Fano defect on its 
velocity V. The parameters are C = 7] = A = e = 1.0, and 
q = 1. The minimum at Vf = 2CqF is quite shallow in this 
case. 

Still we observe a smooth T(V) curve (Fig. |SJ). This 
result implies that the plane- wave modes which make up 
the soliton loose their resonant backscattering features. 
This is not unexpected, as the interaction of many modes 
between themselves can be interpreted as a multichan- 
nel scattering problem (as opposed to many uncoupled 
single-channel scattering processes in the absence of the 
mode- mode interaction). Multichannel scattering prob- 
lems are expected to have less pronounced Fano reso- 
nance features, since the Fano resonance is inherently 
based on keeping phase coherence in order to support 
the destructive interference. As the phase coherence gets 
suppressed by the mode-mode interactions, the necessary 
interference is detuned. 



Resonant transmission of solitons 

The passage of a soliton through a single-site impurity 
has been already studied in many works |2lj . Here we 
consider the nonlinear extension of the two-site-impurity 
system Q, 

i<j> n — C{(t>n-l+(t)n+l) + \\4'n\ 2 (t>n + (Sn,0 + 8n,l)e4'n ■ (23) 

The transmission coefficient for the soliton passing this 
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FIG. 6: A typical dependence of the transmission coefficient 
T for the soliton in the model 1231 on V for rj = A = 1. and 
(a) e = 19 and C = 10 (solid line); (b) e = 3.8 and G = 2 
(dashed line); (c) e = 2.85 and C = 1.5 (dotted line) In this 
figure the normalized velocity V = is used. In all the 
cases, the position of the resonance is fixed at V re s ~ 0.318. 



defect was computed numerically for various values of 
parameters, and found to be in a very good agreement 
with Eq. © , sec FigO Remarkably, neither the position 
of the resonance nor the value of the transmission coeffi- 
cient at the resonance is conspicuously affected when the 
wave- wave interactions become important. This clearly 
indicates that dephasing does not seriously harm the 
resonant-transmission mechanism. 



CONCLUSIONS 



We have shown that solitons may be resonantly trans- 
mitted or backscattered depending on their velocity and 
the type of the scatterer. We explained these effects by 
considering the soliton as a superposition of plane-wave 
modes in the limit when the mode-mode interaction is 
not affecting the scattering process. By tuning parame- 
ters into a regime where the mode-mode interaction be- 
comes essential during the scattering process, we have ob- 
served that the phase- insensitive resonant transmission is 
not affected. However, the Fano-resonant backscattering, 
which relies on the phase coherence in the two-channel 
scattering process, is completely suppressed in this case, 
which is a consequence of the dephasing of individual 
modes due to the interaction between them. 
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